Thermal phase transitions of supersolids in the extended Bose-Hubbard model 
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We investigate numerically the finite-temperature phase diagrams of the extended Bose-Hubbard 
model in a two-dimensional square lattice. In particular, we focus on the melting of supersolid 
phases of two different crystal orderings, stripe and star orders, arising from the competition of the 
nearest- and next-nearest -neighbor interactions in the vicinity of quarter filling. The two crystal 
orders are the result of broken translational symmetry in either one or in both x, and y directions. 
The broken gauge symmetry of the supersolids are found to be restored via a Kosterlitz-Thouless 
transition while the broken translational symmetries are restored via a single second-order phase 
transition, instead of two second-order transitions in the Ising universality class. On the other hand, 
the phase transitions between the star and stripe orders are first order in nature. 

PACS numbers: 75.10.Jm, 05.30.Jp, 75.45.+J, 75.40.Mg 



I. INTRODUCTION 

Supersolids with both diagonal and off-diagonal long- 
range order are observed in various model, either with or 
without hard-core constraint One obvious necessary 
condition for the occurrence of supersolids (SSs) is the 
presence of a small but finite kinetic energy in competi- 
tion with a relatively large repulsive interactions. Previ- 
ous results suggest that the hard-core constraint may de- 
stroy supersolid phases by the formation of domain walls, 
which leads to phase separation instead^ Stable super- 
solid phases in hard-core models are found, however, in 
systems contain frustrated interactions*^— 

The extended Bose-Hubbard model is a typical exam- 
ple that has been proposed to demonstrate a supersold 
phase within the mean-field approximation^, and is con- 
firmed by quantum Monte Carlo simulations As 
long as the next-nearest-neighbor (nnn) interaction V% is 
dominant, hard-core bosons lining up as stripes to reduce 
potential energy near half filling^ (see Fig. [T]). At quarter 
filling, frustrations induced by competing nearest neigh- 
bor (nn) Vi can lead to a stable star order phase that 
characterized by finite structure factors at wave vector 
(tt, 7t) and (w, 0) [(0, 7r)]. The existence of the correspond- 
ing supersolid is also numerically confirmed recently^— 

The melting of stripe supersolids has been investigated 
in Ref. 9 where the finite-temperature phase diagrams 
are determined. At half filling, the stripes melt at a first 
order phase transition, while in the doped system, the 
melting transitions of the solid and supersolid are either 
very weakly first order or of second order. Both the bro- 
ken rotational and translational symmetries of the stripe 
phases are found to be restored at a single continuous 
transition and no intermediate nematic (liquid crystal) 
phase has been observed. Nevertheless, Ref. 9 focuses 
solely on the case with no nn interaction V\ and so only 
the melting of stripe phases are discussed. The thermal 
melting behavior of the star supersolid phases near the 
quarter filling is then largely unknown. Since the star or- 
der is characterized by two structure factors, the two bro- 
ken translation symmetries can, in principle, be restored 



simultaneously or at two different temperatures. The 
latter scenario means the star phase melts into a stripe 
phase before changes to a normal fluid or supcrfluid. It 
is also worth to examine whether the missing nematic 
phase in Ref. 9 can be stabilized by the introducing V\ 
in the system. It is the purpose of this work to address 
these issues and determine the finite-temperature phase 
diagram for the competing nn and nnn interactions. 



II. MODEL 

To start with, the extended Bose-Hubbard model on 
a two-dimensional square lattice with the Hamiltonian is 
defined as 



H b = -tY j {b\b, J +b l b])+V l 



1 ^2 n>i1lj+V2 ^ UiTlj-fl ^2 n i 

with hard-core constraint nj < 1. In the above formula, 
b(tf) is the boson destruction (creation) operator and 
^) "(£"" ) sums over the nearest (next-nearest) neigh- 
boring sites. The energy scale is fixed by t = 1 through- 
out this paper. 

In the case of V% = 0, it is well known that the ground 
state at half filling is a checkerbroad state, which becomes 
a superfluid as the system is doped away from half filling. 
On the other hand, in the case of V\ = 0, the half-filled 
ground state has a stripe order that characterized by the 
finite structure factors at wave vector Q= (n, 0) or Q= 
(0, 7r) (Fig. [T|). More interesting is the coexistence of the 
stripe order and superfluidity, i.e., the stripe supersolid 
phase, in the doped system. Most condensed bosons in 
this exotic states move in the direction parallel to the 
stripes but superflow in the perpendicular direction is 
also observed^ The stripe quantum solid proceeds to the 
supcrfluid via two second-order transitions, in contrast to 
the case of vanishing V% , where transitions from quantum 
solid to superfluid are discontinuous. 

The phase diagram becomes even richer when both Vi 
and V2 are finite. An extra quarter-filled solid (see Fig. 
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FIG. 1: A schematic description of the classical (a) stripe 
order and (b) star order for half and quarter fillings, respec- 
tively. While there is only one broken translational symmetry 
in the stripe phase, two translational symmetries in both di- 
rections are broken in the star phase. On the other hand, the 
x-y rotational symmetry is broken in the stripe phase but not 
in the star phase. 

[1]) and its corresponding supersolid are predicted within 
the mean-field theory-^ and are confirmed by numerically 
simulations.— Quarter-filled solids with two distinct 
crystal structures are possible depending on the compe- 
tition of Vi and V2, while only one structure is stable 
for the supersolid phased This so-called star supersolid 
has been discussed in detail but the finite-temperature 
phase diagram has not been examined so far and there- 
fore serves as the subject of this work. 

One of the most efficient approach to study the bosonic 
models is the standard stochastic series expansion (SSE) 
Monte Carlo method implemented with a directed loop 
algorithm. 10 In SSE, the superfluidity, given by p x (y) = 
(W^/ja) / '4/Jt, is computed by measuring the winding 
number fluctuation within the simulations. Hereafter we 
present the average superfluidity p s = (p x + p y )/2. The 
broken translational symmetries in the x and y directions 
of the striped phases are characterized by the structure 
factors 

S(Q)= y £(n i n j e i <*«)/N*, (2) 

ij 

at the wave vectors Q = (n, 0) and (0,tt), respectively. 
For convenience, we measure the average value of 

S+=[S(7r,0)+S(0,7r)}/2, (3) 

in order to avoid the difficulty to determine the direc- 
tion of broken translational symmetry in each simula- 
tion. Moreover, we also study the rotational symmetry 
breaking in the local-density correlation by measuring 
the order parameter On, 

On = n x , y (n x +i,y - n XtV+ i), (4) 

x,y 

with n XyV the boson number operator at the site (x,y). 
We stress that On signals the local broken rotational 
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FIG. 2: (Color online) Ground-state phase diagram for V\ = 
2.0. It contains a half-filled stripe quantum solid and a 
quarter-filled star quantum solid and their corresponding SSs. 
Finite-temperature phase diagrams scan along the lines A and 
B are presented in Figs. [3] and [H] respectively. The dotted 
line indicates a first-order transition from the star supersolid 
to stripe supersolid. The system size used here is L = 28 with 
T = 1/2L. 

symmetry but not the global one. The star order phase, 
on the other hand, is characterized by the finite values 
of both S(jr, 7r) and S+, which breaks the translational 
symmetry in both directions. 



III. GROUND-STATE PHASE DIAGRAM 

The ground-state phase diagrams V\ vs V2 at quarter 
filling have been presented before^. The star supersolid 
is found to be stable for a wide region of parameters 
V\ and V2 even at commensurate filling n — 1/4. For 
the following discussion, we choose a representative nn 
interaction V\ — 2.0 of which the phase diagram pvs. t is 
shown in Fig. [5] Because of the particle-hole symmetry, 
the symmetric upper half of the diagram is not shown. 
The main lobe in the diagram is the striped quantum 
solid phase at half filling and from which decreasing the 
boson number leads to a striped supersolid phase. 

As mentioned before, the striped supersolid is char- 
acterized by a finite S+ that breaks one translational 
symmetry but preserves the translational symmetry in 
another direction (Fig. [T|). Around quarter filling, the 
£(71", 7r) becomes finite as the system enters the star su- 
persolid phase. The translational symmetry is now bro- 
ken in both x and y directions. At exactly quarter filling 
the superfluidity p s vanishes for V2 > 4.9 (sec Fig. [3]) but 
remains finite otherwise. 

As shown in the Fig. 1 of Ref. 8, the general features of 
the ground-state phase diagram can also be described by 
the mean-field theory. The quarter-filled star quantum 
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FIG. 3: (Color online) Finite-temperature phase diagram at 
quarter filling along the line A in Fig. [2] The black solid 
lines are guides to the eyes. Data points are obtained from 
the finite-size analysis as described in Figs. [4] and [5] In the 
following figures, error bars are smaller the symbol size unless 
otherwise stated. 

solid phase is enclosed by a narrow region of star su- 
persolid phase before enters the stripe supersolid phase. 
However, the mean-field approximation tends to overes- 
timate the extent of the stripe and star supersolid phases 
as presented in the Fig. 8 of Ref. 7. This indicates that 
the neglected quantum fluctuations in the mean-field cal- 
culations destabilize the supersolid phases and reduce the 
region of the supersolid phases in the ground-state phase 
diagram. 

In the following, we present the finite-temperature 
phase diagrams at quarter filling and away from quar- 
ter filling along the lines A and B in Fig. [2j respectively. 

IV. FINITE-TEMPERATURE PHASE 
DIAGRAM AT QUARTER FILLING 

Figure [3J shows the finite-temperature phase diagram 
with fixed boson density of 1/4 (line A in Fig. [5]). In our 
grand-canonical simulation we adjust the chemical po- 
tential p in order to obtain the average particle density 
n equals to quarter filling. At zero temperature, the off- 
diagonal long-range order is depressed when Vi increases, 
while the diagonal long-range order, the star order in this 
case, is enhanced. In the intermediate V%, both orderings 
coexist and gives rise to a star supersolid phase. The 
presence of the supersolid phase is in contrast to n=l/2 
case where the superfluid proceeds to the stripe quantum 
solid via a first-order transition with no intermediate su- 
persolid phase. In Fig. [3J the superfluid changes to the 
star quantum solid via two second-order transitions in- 
stead. For the study of thermal transitions, we carry out 
a finite-size analysis on the order parameters p s , S+, and 
S(n, 7r). 

In a two-dimensional system, the superfluid to normal- 
fluid phase transition can be well described by the 
Kosterlitz-Thouless (KT) transition/ 1 ^ Plotting the su- 




FIG. 4: (Color online) Finite-size scaling of the superfluidity 
p a for p/(Vi + V2) = 0.394. The inset shows the logarithmic 
correction T* = T a -t{1 + l/[2ln{L/L )]}. 

perfluidity p s as a function of T, as shown in Fig. |4] for 
a representative V2, the intersection values of p s (T*) — 
—T* for different sizes L follow the logarithmic correction 
T* = T KT {1 + l/[2ln(L/L )}}. Fitting to this relation 
gives the value of critical temperature Txt f°r infinite 
system. Indeed, the values of T* shown in the inset of 
Fig. 0] follow very well the logarithmic correction and 
confirms the superfluid to normal fluid transition is of 
the KT type. The critical temperature Tkt, as shown 
in Fig. [3J reduces smoothly as V% increases and disap- 
pears for V2 > 4.9. 

Next, we analyze the temperature dependence of the 
structure factors S+ and S(ir,ir). As mentioned before, 
the two broken translational symmetries can be restored 
at two different temperatures or restored at the same 
temperature. To determine the transition temperatures, 
we recall that for second-order transitions, the fourth or- 
der Binder cumulant ratios Ul = 1 — (0 4 )/3(0 2 ) 2 should 
intersect at the same critical point for different sizes 
Figure [5] shows the cumulant ratios of S+ and S (it, tt), re- 
spectively. Indeed, U l meet at the same transition points 
independent of the system size. Particularly, the resulted 
transition points from both cumulant ratios are very close 
to each other and are indistinguishable within the statis- 
tical uncertainties of our data. Our data thus suggest 
a single second-order transition for the melting of the 
star quantum solid to the normal-fluid phase, or a very 
narrow intermediate stripe quantum solid in between the 
phases. 

To further clarify these two scenarios, we suppose there 
are two second-order phase transitions that restores the 
translation symmetries sequentially. We expect both 
transitions are in the three-dimensions Ising universality 
class as each transition breaks a Z2 symmetry. Therefore 
the structure factors should follow the finite size scaling 
law S = L- 2f} ' v F[{T - T c )L l l v \ with (3 = 1/8 and v = 1 
as expected for the Ising class. 13 ' 14 The data of SL 2 ^/ U 
for different L should intersect at the same temperature 
that determines the value of T c . By rescaling the tern- 
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FIG. 5: (Color online) Binder cumulant of S+ and S(n,n) as 
a function of temperature for different system sizes (Vi = 2.0, 
V 2 = 4.6). 



FIG. 6: (Color online) Attempt of the finite-size scaling for 
the structure factor S+ (/J./ (V1 + V2) ~ 0.394) by assuming the 
critical exponents /3 = 1/8 and v = 1 of the Ising universality 
class. At low temperatures, significant derivations from the 
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perature to (T — T^L 1 '", all data will collapse into a 
single line, provided the data are close enough to the 
critical point where the scaling law is valid. As shown in 
Figs. H]and [3 the data of S + and S(n, n) intersect at a 
size-independent critical temperature, respectively. The 
critical temperature of S(n, 7r) is found significantly lower 
than that of the 5+ and therefore suggests an interme- 
diate stripe quantum solid phase. However, the rescaled 
data collapse well only for T > T c but a clear derivation 
from the scaling law is observed for temperatures below 
T c . This failure of data collapse invalidates the assump- 
tion of two transitions of Ising type. Therefore it rules out 
the possibility of a very narrow intermediate stripe quan- 
tum solid phase lies between the star quantum solid and 
normal-fluid phases. The melting of the star quantum 
solid to normal fluid is then likely to be a single second- 
order transition. The system happens not to choose to 
restore the translational symmetry in one particular di- 
rection first and the other direction later, but it restores 
both broken translatonal symmetries simultaneously. We 
will comment on this in the next section. 



V. FINITE-TEMPERATURE PHASE DIAGRAM 
AWAY FROM QUARTER FILLING 

When the system is doped away from the quarter fill- 
ing, the finite-temperature phase diagram is even richer. 
Figure |H] presents the finite-temperature phase diagram, 
determined by the Binder cumulants as discussed previ- 
ously, as a function of chemical potential p with fixed 
Vx = 2 (line B in Fig. [2]). The zero-temperature ground 
states are the superfluid, star supersolid, and stripe su- 
persolid phases, respectively, as the boson number is in- 
creased. The quantum phase transition of the superfluid 
to star supersolid is second order, whereas the transition 
from the star supersolid to the stripe supersolid is first or- 
der. Figure IHl plots the ground-state order parameters as 
a function of /1 from which an abrupt change is observed 
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FIG. 7: (Color online) Finite-size analysis of S(tv, tv) 



at the transition between two supersolids. The histogram 
of On (see the inset) shows the double peaks feature as 
a clear signature for a first-order transition. This dis- 
continuous transition reflects the fact that the star or- 
der is more than simply breaking another translational 
symmetry of the stripe order. It implies the two broken 
translational symmetries in the star order are not inde- 
pendent and therefore cannot be constructed by breaking 
the translational symmetries sequentially via two second- 
order transitions. This also explains the direct melting 
of the star quantum solid observed in Fig. |3 

As the temperature is increased, again, the broken 
gauge symmetry is restored via a KT transition. How- 
ever, in contrast to the Fig. [3] where the critical temper- 
ature T c (p s ) reduced smoothly without being affected by 
the emergence of broken translational symmetries, T c (p s ) 
in Fig. [S] exhibits a sudden drop in the star order phase. 
The T c (p s ) remains rather flat throughout the whole star 
supersolid phase. This sudden suppress of superfluidity is 
thus unlikely due to the increase in particle number alone, 
but may be resulted by the coupling of the broken trans- 
lational symmetries. This is also in contrast to the super- 
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FIG. 8: (Color online) Finite-temperature phase diagram 
(Vi = 2.0, V% = 4.6) as a function of fi along the line B 
in F J - □ 




FIG. 9: (Color online) Order parameters as a function of fi 
(Vi = 2.0, V 2 = 4.6, L = 40, T = 0.02). The parameters On 
are plotted for sizes L — 40 (filled circles) and L = 28 (open 
circles) to demonstrate the finite-size effect. Inset: histogram 
of the parameter Oat at (J,/(Vi + V2) ~ 0.4848 shows a double 
peak feature of the first-order transition between the star and 
the stripe supersolid phases. 



solid melting in most systems^— where the restoration 
of superfluid symmetry and translational symmetry are 
largely independent. This interesting behavior may be 
related to the intrinsic structure of the star supersolid 
and further study is needed to gain more insight of this 
particular phase. 

For larger chemical potential, and particle density as 
well, we found that the stripe supersolid melts into a 
stripe quantum solid and then further into the normal 
fluid through a second-order transition. From the plots 
of Binder cumulants for both S+ and On in Fig. [l"0l 
we observe that the transition temperatures for both 
broken translation and rotational symmetries are within 
our data uncertainties and again suggests a single sec- 
ond phase transition. Therefore no intermediate nematic 
phase is observed. While this finding is consistent to the 
case of vanishing nn interaction V\ , whether larger V\ will 
be able to stabilize the nematic phase is worth to study. 
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FIG. 10: (Color online) Binder cumulant of S+ and On as a 
function of temperature for different system sizes (Vi = 2.0, 
V 2 = 4.6). 

VI. SUMMARY 

We present the finite-temperature phase diagrams of 
the extended Bose-Hubbard model in the vicinity of the 
quarter-filling. The transitions from the super to nor- 
mal phases that restore the gauge symmetry are of the 
KT type in nature. Interestingly, the transition tempera- 
ture exhibits a sudden drop when the star order emerges. 
It differs from the usual supersolid that the superfluid 
T c behaves largely independent of the crystal order of 
the system. Futher exploration of the possible coupling 
between the broken gauge and translational symmetries 
should provide more insight of this behavior. 

For the melting of star order, our result support the 
scenario of a single second-order transition instead of two 
Ising-type second-order transitions. Translational sym- 
metries are restored simultaneously and no intermediate 
phase of stripe order is found. Our result indicates that 
the two order parameters of the star order couple to- 
gether and the state cannot be viewed simply as two 
interwined stripe orders. It will be interesting to in- 
vestigate the effect of a small anisotropy in one dimen- 
sion that lifts the degeneracy of the x, and y symmetry. 
As one translational symmetry is more robust, the criti- 
cal temperatures of the broken translational symmetries 
should be different and the melting of star quantum solid 
might occur as two second-order transitions. The effect of 
anisotropy on the stripe order melting will also be worth 
to study to see if the nematic state can be stabilized near 
the stripe order. 

Acknowledgments 

We thank Y. C. Chen and M. F. Yang for stimulating 
discussions. The numerical computations are performed 
in the Center for High Performance Computing of the 
THU. K.K.N, acknowledges the support by the National 
Center for the theoretical Science and the finanical sup- 



6 



port from the NSC (R.O.C.), Grant No. NSC 97-2112- M-029-003-MY3. 



1 G.G. Batrouni and R.T. Scalettar, Phys. Rev. Lett. 84, 
1599 (2000); F. Hebert, et al.. Phys. Rev. B 65, 014513 
(2001). 

2 P. Sengupta, L. P. Pryadko, F. Met, M. Troyer, and G. 
Schmid, Phys. Rev. Lett. 94, 207202 (2005). 

3 S. Wessel and M. Troyer, Phys. Rev. Lett. 95, 127205 
(2005); D. Heidarian and K. Damle, ibid 95, 127206 
(2005); R.G. Melko et al., ibid 95, 127207 (2005). 

4 Y. C. Chen, R. G. Melko, S. Wessel, and Y. J. Kao, Phys. 
Rev. B 77, 014524 (2008). 

5 K. K. Ng and Y. C. Chen, Phys. Rev. B 77, 052506 (2008). 

6 L. Dang, M. Boninsegni, and L. Pollet, Phys. Rev. B 78, 
132512 (2008) 

7 K. K. Ng and Y. C. Chen, J. Phys.: Condens. Matter 22, 
185601 (2010). 

8 A. van Otterlo, K. H. Wagenblast, R. Baltin, C. Bruder, 
R. Fazio, and G. Schon, Phys. Rev. B 52, 16176 (1995). 

9 G. Schmid and M. Troyer, Phys. Rev. Lett. 93, 067003 



(2004). 

10 A. W. Sandvik, Phys. Rev. B 59, R14157 (1999); ibid 56, 
11678 (1997); O. F. Syljuasen and A. W. Sandvik, Phys. 
Rev. E 66, 046701 (2002). 

11 J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 
(1973). 

12 K. Binder, Z. Physik B 43, 119 (1981). 

13 M. P. A. Fisher, P. B. Weichman, G. Grinstein, D. S. 
Fisher, Phys. Rev. B 40, 546 (1989). 

14 N. Lafiorencie and F. Mila, Phys. Rev. Lett. 99, 027202 
(2007). 

15 K. P. Schmidt, J. Dorier, A. M. Lauchli, and F. Mila, Phys. 
Rev. Lett. 100, 090401 (2008). 

16 K. Yamamoto, S. Todo, S. Miyashita, Phys. Rev. B 79, 
094503 (2009). 

17 M. Boninsegni and N. Prokof'ev, Phys. Rev. Lett. 95, 
237204 (2005). 



